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LIGHT PROPAGATION IN 
GENERALLY COVARIANT ELECTRODYNAMICS 
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Within the framework of generally covariant (pre-metric) electrodynamics, we specify 
a local vacuum spacetime relation between the excitation H = {T>, Ti) and the field 
strength F = (E,B). We study the propagation of electromagnetic waves in such a 
spacetime by Hadamard's method and arrive, with the constitutive tensor density k ~ 
dH/dF, at a Fresnel equation which is algebraic of 4th order in the wave covector. We 
determine how the diff'erent pieces of k, in particular the axion and the skewon pieces, 
affect the propagation of light. file dublinws6.tex, 2002-03-27 
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1. Introduction 

Electromagnetic wave propagation is a very important physical phenomenon in 
classical field theory. In general, the geometrical structure of spacetime as well 
as the intrinsic properties and the motion of material media can affect the liEht 
propagation. In the generally covariant pre-metric approach to clectrodynamicJ3'Bi, 
the axioms of electric charge and of magnetic flux conservation manifest themselves 
in the Maxwell equations for the excitation H = (V, H) and the field strength 
F={E,By. 

dH = J, dF = 0. (1) 

These equations should be supplemented by a constitutive law H = H{F). The 
latter relation contains the crucial information about the underlying physical con- 
tinuum (i.e., spacetime and/or material medium). Mathematically, this constitutive 
law arises either from a suitable phenomenological theory of a medium or from the 
electromagnetic field Lagrangian. It can be a nonlinear or even nonlocal relation 
between the electromagnetic excitation and the field strength. 
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Earlier, we have investigated the propagation of waves in the most general linear 
theory. Here we present some results which hold true for all electrodynamical models 
with an arbitrary local spacetime relation. This is how we call the constitutive law 
if it applies to spacetime ( "the vacuum" ) itself. 

If local coordinates are given, with i, j, ... = 0, 1, 2, 3, we can decompose the 
excitation and field strength 2-forms into their components according to 

= i Hij dx' Adx^, F = i Fy dx' A dx' . (2) 



2. Wave propagation: Presnel equation 

We will study the propagation of a discontinuity of the electromagnetic field fol- 
lowing the lines of Ref.cl, see also Refs.u^B. The surface of discontinuity S is defined 
locally by a function $ such that $ = const on S. Across 5*, the geometric Hadamard 
conditions are satisfied: 

[F„]=0, [d,F,k] ^ q^ fjk, (3) 
[ffy]=0, [^,H^k]^q^h,k■ (4) 

Here [JF] (x) denotes the discontinuity of a function T across S, and Qi :— 9^$ is 
the wave covector. Since the spacetime relation H{F) determines the excitation in 
terms of the field strength, the corresponding tensors fij and /i^, describing the 
jumps of the derivatives of field strength and excitation, are related by 

1 dH- 

h^J ^ l^nrj''^ fkl, with K,/' : = . (5) 

In linear electrodynamics, the components of the constitutive tensor Kij^^ are con- 
stant (or, at least, independent of the electromagnetic field). But in general Kij'^^ 
is a function of the electromagnetic field, the velocity of matter, the temperature, 
and other physical and geometrical variables. Quite remarkably, however, all the 
earlier results obtained for linear electrodynamics remain also valid in the general 
case because whatever nonlinear spacetime relation H[F) may exist, the relation 
between the jumps of the field derivatives, according to (^, is always linear. 
If we use Maxwell's equations (||), then (^) and (||) yield 

e'^'^'q.jhki^O, e'^"' q, fki = . (6) 

Introducing the conventional constitutive matrix 

ijkl _ }_ ijmn kl (n\ 

A — 2 "-mn V ' ) 

and making use of (^, we rewrite the above system as 

X^-'''9,/m=0, e^^^'q^f^i^Q. (8) 
Solving the last equation by fij — qiUj — qjai, we finally reduce to 

r^^'q^qkai^O. (9) 
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This algebraic system has a nontrivial solution for only when the determinant 
of the matrix on the l.h.s. vanishes. The latter gives rise to our covariant Fresnel 
equation 

G'''''{x)m,qkqi^o, (10) 

with the fourth order Fresnel tensor density of weight +1 defined by 

G^'^^ix) ^emnp,er.tnX'""'-('x^'l^^l'=x'^'*" . (11) 

It is totally symmetric, 5'^'^'(x) = 0^'-''^''' (x)j ^^^d thus has 35 independent compo- 
nents. 



3. Spacetime relation: The emergence of the axion and the skewon 
fields 

The quantity Ky'^' or, equivalently, x^-'''^ (x), characterizes the electromagnetic prop- 
erties of the vacuum and is as such of universal importance. The untwisted tensor 
density x*-''^'(a;) of weight -f 1 has 36 independent components. If we decompose it 
into irreducible pieces with respect to the 6-dimensional linear group, then we find 

^^jkl ^ (l)^yfei _^ {2)^i3kl ^ {3.) ^ijkl ^ ^j^j^ 36 = 20 ® 15 © 1 (12) 

independent components, respectively. The irreducible pieces of x are defined as 
follows: 

{2)^tjkl ._ 1 l^^ijkl _ -^khj^ ^ _{2)^klij ^ {3)^i]kl ._ ^[ijkl] ^ 

W^^jkl ^ijkl _ {2)^i]kl _ {i)^jkl ^ (l)^fciij ^ (|-^3^ 

Since no metric is available, we cannot foriiLtraces. The Abelian axion piece 
=: a{pc) e*-''^' has been introduced by Ni Q'Q. A constitutive law (for matter) with 
'•^•'x!^' 7^ (non- vanishing "skewon field") has been discussed by Nieves and 
PalBEJ. It yields P- and CP-violating terms in the field equations. Thus all consti- 
tutive functions in ( p^ can claim respectability from a physical point of view in 
the framework of linear response theory. 

Alternatively, by means of mere contractions, we can decompose 

«:,/' = W«:,/'-f(2)^,>' + (3)^^>' (14) 
= ««:./'+2^[,[^<551 +i«5f,<5j]. (15) 

Here := ku'^^ k :— = Kki^^ , and the traceless piece is ^^'^ := Ki' — j kS^. 
By introducing the pseudotensorial skewon field and the pseudoscalar axion field 

by 
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we ultimately find the basic relations between the corresponding irreducible pieces: 

^ 1^.,™„ (l)^^^H^ (17) 

(2),,ijfc' _ ^ ijmn (2) kl ijm[k n '] klm[i n j] f-lQ\ 

X — 2 "'mn — t Om t >->m j V-'-°>' 

4. Properties of the Fresnel tensor density 

The various irreducible pieces of the constitutive tensor affect the wave propagation 
in different ways. Technically, this can be determined by studying how a certain 
irreducible piece contributes to the Fresnel tensor density (0) and, thereby, 
to the Fresnel equation ( |lo|) which governs the wave covectors. 

We will demonstrate here the general properties of the Fresnel tensor density. 
As a preparatory step, let us take x = + Then, using a compact notation by 
omitting the indices, we have quite generally 

Gix) = + + ^ (Oi + O2 + O3 + Ti + + T3) . (20) 

Here the mixed terms Oa contain one 'i/'-factor and the T^'s two ■i/'-factors. Post- 
poning the symmetrization over i,j,k,l to the very last moment, these terms read 
explicitly as follows: 

(21) 
(22) 
(23) 
(24) 
(25) 
(26) 

Now we are in position to prove the following properties of the Fresnel tensor: 

1) Gi'^^^x) — 0. Proof: By direct substitution of — ae*-''^', we have 

This vanishes upon symmetrization over the indices i,j,k,l. 

2) ^('^'x) — 0- Proof: Using the symmetry properties ( p^ and especially the skew 
symmetry (2)-^4ife' — _ (2)^khj ^ gj^j 

c c , {3)^mnri{3) jpsk{3) lqtu__ (3) rimn (3) skjp (3) tulq 

^mnpq ^rstu X A A — ^mnpq ^rstu AAA 

(3) tuql (3) kspj (3) irmn 

— ^mnpq ^rstu AAA 

„_ (3) mnrl {3) kpsj (3) iqtu 

— ^mnpq ^rstu X A A • 



Oi(0,V,</') 


^mnpq 


^rstu 






^Iqtu 


O2(V',0,'/') 


— ^mnpq 


^rstu 


^mnri 


^jpsk 


^Iqtu 


O3(0») 


^mnpq 


^rstu 


^mnri 




^Iqtu 


Ti{ip,(f>,ip) 


— ^mnpq 


^rstu 


^mnri 


^psk 


^Iqtu 


72 (-0, -0,0) 


— ^mnpq 


^rstu 


^mnri 


jjjjpsk 


^Iqtu 


73(0,-0,-0) 


— ^mnpq 


^rstu 




^jpsk 


^Iqtu 
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Upon the symmetrization over the indices k, /, we then get ^C-^^x) = ^ 5('^-'x) 
or Gi^^^x) = 0. 

3) Gi^^^x + + ^^h) = Gi^^h + ^^h)- Proof-. Let us put = ^^'>x + and 
ip = ^'^^X the formula (pc|). Then we can verify that aU the mixed terms (|2l])-([2^) 
vanish. Indeed, using — ae*^'^', we find: 

Oi = 2a emnpq 0™"" 0'?*" (5^ + SPS^Si + S'^SjSZ) 
= 2a emnpq ('?!>''""^'' f/'''?^'' + f/'™"''' 0''^''^' + 

This is zero when we impose the symmetrization over i,j,k,l. Similarly, we find 

O2 = - 2a erstu {6;5l - S^S^gW^'" 0'"*" = 2a erstu (- 0'"" + r^'" 

and 

O3 = - 2a e„ {Si6l ~ ^«5^)0™"" 0^^^'= = 2a e„ (- 0™"'^ 0^f«'= + 0""'?^ 0^^"=). 

Both expressions vanish when we impose symmetrization over i,j,k,l. The proof 
that all T's are equal to zero reduces to the above formulas in which we simply 



need to replace one of the factors by the Levi-Civita e. Consequently, ( |20[ ) yields 

^(0 + yj) = a(0) + GW = ^(0), since GW = Gi^^'h) = 0. 

4) ^/('^^X + ^^^x) = 0. Proof: Take = '^^x and = '^^X- Then the proof reduces 
to the above case. [Formally, we can simply consider a special case of the formula 
^('■^^X + ^^""x + ^^^x) = 0('"^^X + ^^^x) by putting '-^^x = and using the property 



5. Discussion and conclusion 

The properties of the Fresnel tensor simplify greatly the final analysis of the light 
propagation in the general electrodynamical theory. 

As we saw, neither the axion piece '^^^x nor the skewon piece '^^^x (alone or 
together) can provide electromagnetic wave propagation. 

The presence of the principal part ^^^x of the constitutive tensor is indispensable 
for the existence of nontrivial electromagnetic waves and ultimately for the existence 
of the light cone structure on spacetime. In any case, the Fresnel tensor density reads 

^.jfei(^) ^^.j-fei ((1)^^(2)^)^ (27) 

Furthermore, in general 

g^,M((i)^^(2)^)_^^^,M((i)^)^ (28) 

which means that the skewon field does influence the Fresnel equation and thus, 
eventually, the light cone structure. An example of this general result can be found 
in the asymmetric constitutive tensor studied by Nieves and Pal BIIj. 
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Actually, we can use (|^) to find more exactly the contribution of the skewon 
to the Fresnel tensor. One can straightforwardly see that 

Oi((i)x, ^'h, ^'h) = 0, 0,eh, ^'h, ^'^X) + O^eh, ^'h, ^'^X) = 0, (29) 

T2((^)x,(^)x,(^^x)-r3(Wx,(^)x,(^)x). (30) 

The proof is analogous to the demonstration of the property ^/(^^•'x) = ^iid is 
based directly on the skew symmetry (2);^*^'=' — — (2)T^fc'u_ As a result, (EG) yields 



Substituting here the explicit form of the skewon part (18), we finally find 

gy'='(x) = g''''\^^'>x + ^^^x) = ^''■'■'''(^'^x) + (^'x"^*'"'^'^™^™'^ ■ (31) 

Our basic results on the light propagation are then that ^^•'x ^ 0, otherwise 
there is no orderly wave propagation, the axion piece a is left arbitrary, it doesn't 
influence light propagation locally, whereas the skewon piece Si^ =9i-' makes itself 
felt by means of equation (|3l|). The detailed study of how exactly these 15 functions 

affect the light cone is left for future work. 
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